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I. INTRODUCTION 



An interesting generalization of the Boltzmann-Gibbs entropy form has been recently proposed by Tsallis Q , for 
a recent review see for a full bibliography see Q. This new entropy that we discuss below, posseses the 

usual properties of positivity, concavity and irreversibility and generalize the additivity in a non-extensive way. In 
the later years, a lot of work has been done in order to show either that many relationships involving energy and 
entropy in the usual scheme holds in the new one or in order to search for suitable generalizations for them. We 
should mention, among others, the work of Curado and Tsallis ||, that shows that the whole mathematical structure 
of thermodynamics -Legendre transform based- becomes invariant under a change from extensive to non-extensive 
CD 1 statistics (NES). This, in fact, is a property recently proved to be valid for any entropic functional form making 
then important to have a partition function within the formalism. Boltzmann-Gibbs thermostatics constitutes very 
complete and powerful techniques in situations whenever thermodynamics extensivity holds. But this formalism 
provide divergent partition, energy and entropy functions whenever the physical system includes long range forces 
or long memory effects. In this regard, the connection with quantum mechanics and with information theory has 
been stablished |tJ and further applications to self-gravitating objects [|| and astrophysics Q, to Levy flights |lO[ , 
turbulence of fluids |IlJ and to the problem of solar neutrinos jl^] have also been studied. In general, systems which 
present a fractal like or unconventional structure in their space-time description or in their phase space evolution could 
(— > ■ develop hardly tractable mathematical problems within the standard formalism || . Very recently, Grigera pl| derived 
a molecular dynamic test for systems of long range interacting particles with potentials of the form A/r 12 — B/r a . 
He found that the potential energy per particle do not converge for some cases of a and sizes of the system, whereas 
non-extensive quantities do. Also see Refs. HH- 

Non-extensive thermostatistics was also used to study the cosmic blackbody radiation Jl6|,|l7j . This was the setting 
in which powerful bounds upon the non-extensive parameter were obtained. We shall review some of these results 
in our final section. A new analysis of NES blackbody cosmological radiation showed that, independently of the 
degree of non-extensivity, the temperature of the universe varies as the inverse of the scale factor of a flat Friedmann- 
Robertson- Walker metric []l8| . This result suggested cosmology of the early universe as a very profitable arena for 
searching observational consequences of NES. As far as we know, the first of such studies were presented in Jl9[] , where 
a stringent bound on non-extensivity, of comparable order of magnitude as those obtained via the cosmic blackbody 
radiation, was derived. There, it was analized how a first order variation in the non-extensive parameter would affect 
the amount of primordial Helium production, a quantity which can be compared with observational results, via a 
rather easy and straightforward approximate technique. 

We now turn to the issue of early universe cosmology and particularly, to the computation of primordial neutron to 
baryon ratio in the framework of NES. We shall also analyze cosmological situations of great importance, such as pair 
annihilation, temperature relationships, and comoving quantities. We thus refine the previous study and enhance it 
by a sharpening of the model. Although we shall present a very brief review of the main formulae of NES; in first 
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approach, this work should be understood as the search of deviations from the standard model of the early universe 
provided by a different statistical description. 

What we want to use is a direct analytical computation of Helium production, following the work of Bernstein, 
Brown and Feimberg for the standard model, to allow for an approximate study. It is a matter of fact, that 
even when nucleosynthesis processes may be as complicated as one can afford, some simple assumptions makes results 
agree, within a typical few percent, with those generated by a numerical code, thus allowing for the isolation of the 
essential physics. That is exactly what it is needed here in order to study how, for instance, a different formula for 
the number of particles influence the evolution of the early universe and how can this be used to extract testeable 
predictions. Nucleosynthesis has given a strong basis to do such a thing when alternatives theories of gravitation are 
concerned, see for instance ]2l|-p3|], and also for particle physics [^4]-^6] and we hope to show, the same is applicable 
for alternative statistics. 

The rest of the work is organized as follows. Firstly, some useful concepts and formulae of NES are presented. 
Afterwards, the analytical computation of the neutron to baryon ratio is done, and finally we study the necessary 
approach to handle with neutron decay corrections. In the final section, we can use our previous results to estimate 
the mass fraction of primordial Helium produced in the context of NES. Although we shall not have at that moment 
every ingredient needed to do such a thing with great precision, the bound which we obtain hardly constrain the 
degree of non-extensivity in the very early universe and is already comparable with previous results. Finally, we state 
our conclusions and mention other recently obtained related bounds. 



II. NONEXTENSIVE STATISTICS 



Let us first recall some formulae of nonextensive statistic (NES), in Tsallis's approach. The formalism starts by 
postulating |l| 

Postulate 1.- The entropy of a system that can be found with probability pi in any ofW different microstates i is 
given by 



= kjz TT >> " Pi\ = k — l-y>? I > (!) 



w / w \ 

with q a real parameter. We have a different statistics for every possible q— value. In (Q) we have used, of course, that 

5> = 1 - (2) 



In general, 



S q = k'-^f. (3) 



Postulate 2.- An experimental measurement of an observable A, whose expectation value in microstate i is ai, 
yields the q- expectation value (generalized expectation value (GEV)) 



< A > q = Y / p' 1 i a i = Trp«A, (4) 



for the observable A. 

These two statements have the rank of axioms. As such, their validity is to be decided exclusively by the conclusions 
to which they lead, and ultimately by comparison with experiment. The entropy given in (fil) is non-negative and 
reproduces the Boltzmann-Gibbs one (Si = k J2i Pi m Pi) in the limit q — ► 1. In fact, note that (P can be obtained from 
the Boltzmann-Gibbs expression by writing ln(x) = lim g ^i(a; 9_1 — l)/(<z — 1) and then dropping the limit operator. 

It satisfies the pseudoadditivity property which states that if A and B are two independent systems (p[f +B ^ = pfpf) 
then, 

S q (A + B) _ S q (A) S q {B) S q (A) S q (B) 

T, — Z 1 T, r(l-q) — r r — , (oj 
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which shows that (1 — q) is a measure of the degree of non-extensivity of the system. The optimization of S q given 
by @ together with the constraints Trp = 1 and Trp q H — U q < oo yield to the canonical ensemble equilibrium 
distribution 



1 r 



1 - (1 - q)0H 



(6) 



and to the generalized partition function, 



Tr 



i-{i- q )m 



(7) 



Here, as usual, (3 = 1/kT and H is the hamiltonian of the system. We shall focus in the (3(q — 1) — > limit, in which 
a first order expansion allows analytical computations. The expression of a general mean value of an operator was 
computed in this limit by Tsallis, Sa Barreto and Loh pjjfl. When applied to particle number operators, the result is: 



< h >„=< n >bg Z 



9-1 

BG 



l+(l-q)x 



< tr >bg 

< n >bg 



. » 2 ^ < n" > BG 
< n >BG -^r— 

< n >bg 



(8) 



where x stands for e/kT (e is the energy of a simple particle) and the symbol BG means to be computed in Boltzmann- 
Gibbs statistics. With the standard values of < h 2 >bg an d < n 3 >bg, f° r fcrmions and bosons, we obtain the main 
corrections (the second terms in equation (8)), Cb OS ons,fermions, as: 



C, 



bosons 



1 + cxp (- 



x 1 — |— 3 exp (—x) 



exp (x) — 1 \ 1 — exp (— x) 2(1— exp (— 



(9) 



C 



f ermions 



exp(x) 



x ( 1 1 + exp (2x) + 2 exp (x) 

1 X \ + 2 V 1 + exp (x) 



(10) 



When x is such that exp (x) ^> 1, then, Cb osons 

Cbosons — Cf erm ions — ( 2* /2) exp ( — ,x). 



f ermions 



(1 + cxp (a;)) 2 
x(l — x/2) exp (— x); while for x big enough we get 



III. NEUTRON-PROTON RATIO IN AN EXPANDING BACKGROUND 



We turn now our attention to the description of the model for the evolution of the neutron abundance as the 
universe evolves. As stated in the introduction, we shall follow the leading ideas of Bernstein, Brown and Feimberg 
pfH . We shall denote by X pn (T(t)) the rate for the weak processes to convert protons into neutrons and by X np (T(t)) 
the rate for the reverse ones. X(T(t)) will be, as usual, the number of neutrons to the total number of baryons. For 
it, a kinetic equation may be built, 



dX(t) 
dt 



X pn (T)(l-X(t))-X np {T)X(t), 



The solution to the previous equation reads: 



X(T)= / dt'l{t,t')X pn (t')+X(t )I(t,t ), 



with the integrating factor given by, 



where, 



I(t,t') = exp (- J 



diA{i) 



A(i) = X p „{t) +X„ p (t). 



(11) 



(12) 



(13) 



(14) 



Note that the form of this solution does not depend on the statistical mechanics used. However, changing the statistical 
mechanics entails modifications in the rates, and this will make the solution differs from the standard. Expression 
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( |l2| ) is simplified by setting the initial value to = 0. We expect that X np and X pn be very large at early times and 
high temperatures. Then, the integration factor J(t,io) will be very small for times order t ~ l/A(io) and therefore 
the term X(to)I (i, to) can be omitted. Besides, we also expect that due to the large interaction rates, the integral in 
the first term will be not sensitive to the change of to by 0. With this choice, we now have, 



t 

>TU 4.1 



X(t) = / dt'l(t,t')X pn . (15) 



Finally, we may note that 



and integrate by parts to obtain 



X(!) = w-/ V(u '^(w)- (17) 



A. Rate Formulae 

To explicitly compute expression (|l7|), we need to know the functional form for the rates. Let us first consider the 
rate X np (t) . This is the sum of the rates of three individual processes, 



which are given by 



Xnp ^u+n^p+e+ F <^e++n^p+£' "F ^n^p+e~ +0 (1^) 
/>oo 

K+n^p+e- = A j dp v p 2 u p e E e (\- < h e >) < h v > (19) 
Jo 

poo 

K++n^p+i> = A j dp e p 2 e p u E u (l- < n v >) < h e > (20) 



rpo 

K^p+e-+u = A / dp e p\p v E v (\- < h v >)(1- < n e >), (21) 
Jo 

where A is an overall effective constant -fixed by the experimental value of X n ^ p+e - +£>-, p v ,e stands for the magnitudes 
of the neutrino and electron momentum and F?„ >e for their energies. The recoil energy of the nucleons may be neglected, 
and this enable us to write the energy conservation equation E u + m n = E e + m p for ( |l9]) and E u F m p = E e + m n for 
(pp|). These equations must be used in order to explicitly compute the integrals. In (]2l|), E v = Am — E e > 0, -with 
Am = m n — m p — 1.29MeV- and this gives the upper limit in the integration range. Finally, n„ je are distributions 
functions and (1 — n v _ e ) are blocking Pauli factors. In the usual scheme, n„. e are given by Fermi distributions, whereas 
in our case, by the < n v . e > q — n^ e (q) of the previous section. In general, the electron and neutrino temperatures, T e 
and T v , may differ because at the end of the freezing out period, electrons and positrons annihilate, heating only the 
photons and mantaining with them thermal equilibrium. The magnitude of this difference will be discussed below. As 
we shall see, this difference will be of a few percent during the period of consideration and our first approximations, 
following Bernstein et. a!., is to set all temperatures equal, T = T e = T v = T 7 . When working in the Boltzmann- 
Gibbs framework this unique assumption ensure that the rates for reverse reactions, such as e~ + p — ► n + v, obey 
the principle of detailed balance (PDB), that we analize below in this general setting. 

B. Principle of Detailed Balance 

Let us consider, in the Boltzmann-Gibbs (BG) statistics and as an example, the reverse rate for (fil)]), that is, 
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p oc 

X e - +P ^ n+1/ = A dp e p 2 e p v E v (l- <h v >) < n e > . (22) 
Jo 

Using energy conservation and the fact that p e dp e = E e dE e we have, neglecting the recoil energy of the nucleon, 
dE e = dE v = dp Vl and thus: 



/>oo 

K-+ P ~>n+v = A dp v p\p e E e (\- < fi u >) <n e > 
Jo 



(23) 



As < n„ e >= n u e are Fermi-Dirac distributions, it is easily seen that (1 — n v ) = 7i„ exp (_E„/T) and that n e = 
(1 — n e ) exp {—E e /T). Replacing this into ( p3| ) and recalling that Am = E e — E v , we get 

\ e - +p ^ n+v = exp (-Am/T)\ n+u ^ e - +p (24) 

which is the expression of the PDB. 

When NES is concerned, and correspondingly change from the Fermi-Dirac distributions to the n(q) ones, the PDB 
is no longer valid. In this framework, it may be seen that 

exp(-Am/T)(l - n e {q))n v {q) = (1 - n v {q))n e {q) + (q- 1)0, (25) 

where the correction is 

/ / Am\ ( E e 

4> = (1 - n ei BG)C v I exp I — — I +exp I — — 

n v ,BGC e ^exp (-^p^J + exp (^f^) ( 26 ) 

and the C factors are as in (|l0|). In order to get a simplified picture, and following the BG case, we further approximate 
the model -exactly as was made in the standard case- by assuming that during the period of freezing, the temperature 
T is low in comparison of the typical energies E that contributes in the integrals for the rates. Hence, we may 
replace the Fermi-Dirac distributions by Boltzmann weights (i.e. n„ je ~ exp (— E/T)) and consistently neglect the 
the Pauli blockings (i.e. (1 — n u-e (q)) ~ 1). Even in the case of nonextensive statistics, Pauli Blockings corrections 
are 1 — n q = 1 — tibg(1 — x 2 exp (— x)(q — l)/2), which are neglectable for x S> 1 and a first order deviation from 
q = 1. In addition, doing this allows for the analytical study that follows and which is the objective of the work. In 
this i» 1 regime, the one which we are going to use thorough, we obtain C e ~ C v exp (— Am/T), and we recover for 
the reverse rates the expression given in (p4|) of the PDB. 



C. Rates Computation 

Before we finally compute the integrals for the rates we neglect the electron mass in equations ( pj| ) and ( pp| ) in 
comparison to the energies E v<e . When this is made, the two rates become identical, exactly as they were in BG 
statistics. Then we obtain: 

K+n^p+e- = K+nXt- + ( 480T5 + 2 X 84T ^ + 18T 3 Q 2 ) (1 - q)A (27) 

where, 

K+nX+e- = ( 4!T2 + 2 X 3!T< 3 + 2!< 2 2 ) AT3 - ( 28 ) 

Concerning the rate of free decay of neutrons, we note that NES does not present any change within the context of 
previous approximations. That is due to the disappearence of all distributions functions from (21). In that equation 
we cannot neglect the mass of the electron and the standard result holds, 

- = K^ P +e-+u = 0.0157^1Am 5 . (29) 
This enable to eliminate A in favor of the measured quantity r, 

A = --Am 5 , a = 255. (30) 
t 4 
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All along in this section, we shall neglect the free neutron decay rate when computing the total X np . Thus, we have 

Anp — 2A„+ n _>p+ e - . (31) 
Finally, by using the variable y — Am/T, we get 

A «p = t( 12 + 6i j + y 2 ) + (!-?) ( 240 + M y + Q y 2 )] > ( 32 ) 

for the neutron-proton rate and 

X pn = cxp {-y)\ np , (33) 

for the reverse reaction. 

D. Evolution of the neutron abundance 

In equation (|l7|), we change now variables to the scaled temperature y and obtain, 



i = Apn(y) [ V 1.JTI.. ../^ d f \n(y') 

The integrating factor now becomes, 



/>(!) a( « 



J(y,y')=exp - / dy — A(y) . (35) 



To evaluate the jacobian dt/dy, we need to recall that the scale factor of the universe, i?, in a Friedmann-Robertson- 
Walker metric, behaves as R ~ 1/T, independently of the statistics p^| . Therefore, T/T = —R/R, and the right 
hand side is given by Einstein equations: 



R /8ttG \ 

R = {— p J 



(36) 



Here p is the energy density for relativistic particles in NES. This may be easily computed, if the distributions functions 
are known, by 

/>oo 

P=^J o E 3 n(q)dE, (37) 
where g is the degeneracy factor. The result is 

7T 2 : 1 

P = Pbosons + Pferrnions = T^gT 4 + — ~ (40.02 56 + 34.70 g f) T A (q - 1), (38) 

where g = gb + 7/8gj. At high enough temperatures the energy density of the universe is essentially dominated by 
e~ , e + , v and 7's. Interactions among these particles keep all them nearly the same temperature. Accordingly, we 
set, <?b = 2, g/ = 2 + 2 + 2x3 = 10 and g = 43/4. With these values we get, 

p=^ 5 T 4 + ^21.63T 4 (g-l), (39) 

With this in mind, and using T/T = —a/a, we obtain 

'4ir 3 Gg\ (AG 



dT _ 
~dt ~~ 

which yields, to first order in (q — 1) to, 



45 / V 3tt 



21.63(9- 1) 



T 3 (40) 
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We shall call the constants, 



_ / 45 Y 1 



and 

45 \ 21.63 



= b 



67r 4 / g 



With all above, the integrating factor becomes, 

I(y,y')=exp(K(y)-K(y% 

where, 



K(y) = -b dy 



12 6 1 \ , , u /240 84 9 \ , 

y 4 y A y 2 J \ y 4 y A y 2 J 

(1 + exp(-y)). 

This integrates to give, 



/80 


42 


9 


/80 2 7\ 


W~ 


h — ' 

y 2 


1 h 

y 


\y 3 y 2 yy 



where i?£ stands for the exponential integral [J Introducing, 

^pn(y) _ 1 



Xeq A(y) l + cxp(y)' 
the neutron abundance ratio reads: 

X(y) = X eq + f dy' exp (y')X eq (y') 2 exp [K{y) - K(y% 

J O 

The exact previous result may be expanded again to first order in (q — 1). This gives, 

X(y) = X st 

andard + (1 — q) exp (bA(y))(bB(y) + bcA(y)) x 
V dy'exp(y')X eq (y') 2 exp(bA(y')) 

-(1-q) exp (bA(y)) I dy' exp (y')X eq (y') 2 exp (&A(y')(W?(y') + bcA(y')), 

J O 

where, 

-^standard = X eq (y) + exp (bA(y)) / dy' exp (y')X eq (y') 2 exp (i)vl(y')), 
and the functions A and _B are: 



1 The exponential integral is defined by Ei(n, x) = Gxp t ^„ xt ^ dt for Re(x) > 0. 
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My) = [ 4 + 4 + - + ( 4 + 4 ) ex P i-v) ) (5 1 ) 



o i 9 1 I ^ 9 

t y y \y y 

B(y) = ( ^ + - 2 + - + ( ™ + 4 + 7 ~) (~y) - 7^(1, j/)) . (52) 

We can numerically compute the integrals that appears here. Using explicit values for the constants b = 0.252 and 
c = 0.00244 with the mean life of the neutron given by t = 889.8s we get the curves of Fig. 1. There, it is shown the 
standard value of the neutron to baryon ratio and the correction due to NES. In Fig. 2. are shown different explicit 
computations for a range of the q parameter. Finally, it may be seen that X(y) asymptotes to 

X(y = oo) = 0.15 + (q- 1)1.15, (53) 

for very low temperatures. 



E. Possible dependence on the number of neutrino types 

We briefly assess here the possible dependence on the number of neutrino types. To do this we note that such a 
variation would introduce changes in the coefficients b and c quoted above. For three neutrino families, g — 43/4 and 
the addition of one more family change it in an amount Sg — 7/4. Then, the standard case coefficient b is affected as, 

and the non-extensive related as, 

c 2 g 

Introducing these new values in ( p6[ ) and making all numerical computations again, it is possible to obtain the 
corrections due to a different number of neutrino families. In the sake of conciness, we shall explicitly skip these 
computations. Instead we note that, while in the standard case, the addition of one more neutrino family would 
involve an increment of the energy density of the universe and a speed up in its expansion, in the general case, the 
correction will depend on the sign of the factor (q — 1). 



IV. TOWARDS COMPUTING NEUTRON DECAY CORRECTIONS 

We have already solved the evolution of the neutron abundance neglecting the neutron decay. As in the standard 
formalism, it is useful to change notation by using an overbar for the result just obtained, X(y) — > X(y). Including 
the effects of the free neutron decay in the rate equation we would get, 

X(t) = exp(-t/T)X(t), (56) 

because X does not vary much during the period in which neutrons decay. In the capture time t = t c , when the 
temperature drops somewhat below the deuteron binding energy [ejj — 2.62MeV), the neutrons are captured in 
deuterons. Then, deuterons collide and almost all neutrons present at t = t c are converted into 4 He. Inserting this 
time in (|56j ) and using the asymptotic value of X(y) just derived would yields, thus, half of the mass fraction of helium 
produced in the early universe. In order to obtain a precise value of t = t c or t Ctq for the nonextensive result, we have 
to analize the set of the basic reactions: 

n + p-<=^> L> + 7 (57) 
D + D p + T (58) 

T + D <=^> 4 He + n. (59) 

Doing this in the standard case [^0| yields a value of t c such as the correction exp (— t c /r) ~ 0.8. Although the analysis 
of (^,^,^) in a NES framework is beyond the scope of the present work, we shall establish here the basis on which 
this could be done in the future. For this, we shall study the process of annihilation of pairs e~, e + . 
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A. Conserved comoving quantities 



We know that, when the temperature of the universe is low compared with the electron mass, electrons and positrons 
annihilate heating the photons and, morover, changing the effective quantity of degrees of freedom. Then, one of the 
assumptions of the previous section, i.e. T = T e = T v = T 7 , is no longer valid and T 7 will be different from the 
decoupled neutrino temperature. We shall need the connection between these temperatures. To compute it, we begin 
by noting that, at these times, the number of neutrinos in a comoving volume is fixed, independently of the statistics. 
Effectively, after the moment in which the collision rate for neutrinos is lower than the expansion of the universe, 
R/R, they cannot be neither created nor destroyed. As the number of neutrinos is given by, 
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the relationship, 



or, equivalcntly, 



E < " ><?= 4^2 C(3)ft,T 



n„(t)R(t) 3 = const. 



T v (t)R(t) = const. 



( 9 -l)^5.50T 3 



continues holding. 

/From the full set of Einstein field equations, it may be established that 

d(p(q)R 3 ) = -p(q)d(R 3 ), 

this may also be written as, 



(60) 



(61) 



(62) 



(63) 



dp(q) 
dT ' 



But we may compute p{q) , as was done with the energy density, obtaining p{q) 
this and (|64[), we get 



d_ 

dT 



p{q) +p(q) 

T 



R< 



o. 



(64) 

p{q)/3 for the relativistic gas. Using 



(65) 



as can be cheked by direct differentiation. This means, that in general case, (p(q) +p(q))/T is a conserved quantity 
in a comoving volume. When q = 1 the previous expression stands for the entropy of the system. 



B. Pair annihilation 



We may now consider that pair annihilation occurs at a temperature T, and then indicate with symbols — and 
quantities before and after T. From all above, we have 



p{q)- +p(q)- = p(q)++p(q)+ 



(66) 



Here p(q)- t + is as in ( p8| ) with the corresponding values of <?-,+ and gb.j\- t +. Due to the removal of pairs, we have 
that g+ < <7_. We may now expand (p6|), to first order in (q— 1), obtaning the relationship, 

Tl = T 3 _ (j^j [l + (q- 1)6} , (67) 

where, 

30 30 
6 = (40.02 56 ,_ + 34.70 ff/ ,_)— (40.02 ffb ,+ + 34.70$/,+) 5-5— . (68) 

Z7T <7_ Z7T <7+ 

Until the moment of e - , e + annihilation, the gas composed by e~, e + and 7 follows a law of identical form that (|6^), 
as such temperatures, being T e = T 7 . But after the annihilation, where only 7's remain, we have 

T 7 - T v (^j 3 [1 + (1 - g )0.013] , (69) 
where the corresponding values g- = g(e~ , e + , 7) = 11/2 and g+ = g(j) = 2 were used. 
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C. Neutron capture time 



To finally compute the capture time t c , we use that 

1 dT v /8ttG 



T v dt 



\ 3 



from where 



t = 

follows. Here p(q) adquires its low energy form 

■2 



PO 



aJ-^Tt + (q - 1)34.707^ 



Tv T' v \8nG P (q) 



5 7 ^T 7 4 + (g-l)40.02T^ 



where </„ = 3 and g 1 = 2. Recalling the relationship Q6S|) we get, to first order in {q — 1 



(70) 



(71) 



(72) 



Po=g e ff—T*+(q-l)TtUM 



(73) 



where g eff = ,g„ + ,g 7 (ll/4) 4 / 3 . 

In the neighbourhood of the time of the neutron capture , an approximate evaluation of J7l] ) is in order. Using the 
expression for po just derived, we get 



/ 45 



\l6^g eff Gj T* 



r (l + (1-5)1.74)+ t , 



(74) 



where to is an additional integration constant and the corrective number 1.74 arise from the first order evaluation of 
p{q)z. Recalling then the relationship between the temperatures, we may express 



t = 



45 



16ir 3 g ef fG 



— (1 + (1-3)1.76)+ to, 



(75) 



We may neglect to in a first approximation or obtain it from a perturbation analysis of ( |7l"| ) , [ p0| . 

Now, a comprehensive analysis of the reactions (|57],|58|,^9]) would yield a value for T 7 , for which almost all nucleons 
are converted into helium . Using then (f75|), the evaluation of t c follows and thus, the mass fraction of primordial 
helium. In passing through NES, we have to note that a fundamental tool in this analysis, i.e. Saha formulae, is 
modified in a non-trivial way. This entails modifications in the abundance fractions that may in turn change their 
ratios, which are the magnitudes on which the set depends. So, although phcnomenological fits for the rates may be 
considered, we acknowledge here that further research must be done to analytically solve for the reactions, in the case 
this be ultimately possible. This analysis is far out from the present work and we shall not discuss this problem in 
deeper detail. 



V. CONCLUSIONS AND COMPARISON WITH OTHERS RELATED RESULTS 

We have presented a comprehensive study of the primordial neutron-baryon ratio in the context of non-extensive 
statistics. This work, thus, extends previous advances on the helium nucleosynthesis production problem [ fl9f and is 
the result of a detailed analysis of the weak interaction rate and the energy densities corrections at first order in the 
non-extensive parameter (q — 1). We have also studied some important cosmological scenarios such as the possible 
existence of conserved comoving quantities, temperatures relationships and pair annihilation, together with pointing 
out where further research is needed. These other results will be used elsewhere to put new stringent bounds on the 
non-extensivity of the universe by using high precision satellite cosmological measurements of the cosmic microwave 
background and new galaxy surveys. 
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Even without getting a precise analytical value for the Helium synthesis, it is possible to obtain a first insight in 
the problem by considering it as exactly twice of the neutron-baryon asymptotic value, weighted by the free neutron 
decay: 

ytheoretical _ ^ ^ ( 7 g) 

where X(t) is given by equation (56). But as mentioned above, in order to put an exact bound we would need the 
value of t c q . However, even considering that T 7 is the same as in the standard formalism (T 7 ~ 0.086MeV) we could 
use ( [75| ) to get the correction in t c . Using afterwards the first order expansion of ([36]) (which contains equation (53)) 
we would get the stringent bound: 

|1 - q\ < 2.6 x 1(T 4 . (77) 

In order to get this bound, we have imposed that the non-extensive corrections (coming from equations (75) and 
(53) and the theoretical value of Y p as explained above) be less than the observational error in the primordial 
Helium production, which is doubled so as to neglect the difference between the central values of the analytical and 
observational magnitudes ^4|, i.e.: 

yobservaUonal = q 33 ± q q 2 ( 7 g) 

This abundance is a very safe estimate of the primordial Helium production, which in fact is a very complex observable 
magnitude, see for instance the recent review by Schramm and Turner Q for details. Besides, it encodes other possible 
uncertainties such as the value of the cosmological constant or the exact number of neutrino species. 

The bound ( |77j ) is comparable with that already obtained in Ref. ]l9| . Although it should be recalled that such 
derivation was different from this and the method followed skip the uncertainties of the neutron capture time here 
analyzed. In addition, the weak interaction rates were only approximately treated in | jl9| while here they were studied 
in full detail. 

In a recent letter J28| , it was suggested that the universe as a whole could behave with a non-extensive parameter 
roughly equal to 0.54. Such a value seems to be too high as to resist early universe tests as the one presented here. 

Therefore, this bound, obtained in the simplest situation so as to admit an analytical solution, pertains to those 
imposed in the context of early cosmological models, and from this point of view, enlarge the problems that non- 
extensivity would have as a correct description of the physical universe thorough all eras of cosmic evolution. Very 
recently, other works have also been devoted to study possible deviation from the Boltzmann-Gibbs statistics in 
cosmological and other situations. Among them, and apart from Ref. |lj| already commented, we should note the 
following ones: 

• Plastino, Plastino and Vucetich Jl7j studied the value of the Stefan-Boltzmann constant in the new framework. 
A comparison with the experimental value leads to \q — 1| < 0.67 x 10~ 4 . They also used data from the 
FIRAS experiment on the COBE satellite concening the microwave background radiation, from there they 
obtain \q- 1| < 5.3 x 10~ 4 . 

• Tsallis, Sa Barreto and Loh |l6| also generalized the Planck radiation law and used satellite data from COBE 
to constrain the non-extensive parameter. They found | <gf — 1 1 < 3.6 x 10~ 5 . 

• Tirnakli, Buyiikkilic, and Demirhan [29|j3l[| reworked the generalized quantal distributions functions and also 
obtain a similar bound: \q — 1| < 0.41 x 10~ 4 . 

These works, among others, suggest that non-extensivity is highly constrained concerning the physics of the early 
universe. 
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FIG. 1. Behaviour of the primordial neutron to baryon ratio in standard Boltzmann-Gibbs statistics (curve (1)) and correction 
due to non-extensivity weighted by (q — 1), (curve (2)). The complete behaviour of the neutron-baryon ratio in NES is obtained 
making the sum of curve (1) plus (q — l)x curve (2). 
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FIG. 2. Neutron to baryon ratio relationship for different non-extensive parameter q. 
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